HU-EP-09/51 



43 



43 



- 1—1 

X 



c3 



I. INTRODUCTION 
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chemical potentials 
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The properties of the two-flavored Gross-Neveu model with nonzero current quark mass are inves- 
tigated in the (l+l)-dimensional spacetime at finite isospin /i/ as well as quark number /i chemical 
potentials and zero temperature. The consideration is performed in the limit N c — > oo, i.e. in the 
case with an infinite number of colored quarks. In the plane of parameters fii, [i a rather rich phase 
structure is found, which contains phases with and without pion condensation. We have found a 
great variety of one-quark excitations of these phases, including gapless and gapped quasiparticles. 
Moreover, the mesonic mass spectrum in each phase is also investigated. 
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During the last decade great attention was paid to the investigation of the QCD phase diagram in terms of baryonic 
as well as isotopic (isospin) chemical potentials. First of all, this interest is motivated by experiments on heavy-ion 
collisions, where we have to deal with dense baryonic matter which has an evident isospin asymmetry, i.e. different 
neutron and proton contents of initial ions. Moreover, the dense hadronic/quark matter inside compact stars is also 
isotopically asymmetric. Generally speaking, one of the important QCD applications is just to describe dense and hot 
baryonic matter. However, in the above mentioned realistic situations the density is rather small, and weak coupling 
QCD analysis is not applicable. So, different nonperturbative methods or effective theories such as chiral effective 
Lagrangians and especially Nambu - Jona-Lasinio (NJL) type models [l[ are usually employed for the consideration 
of the properties of dense and hot baryonic matter under heavy-ion experimental and/or compact star conditions, i.e. 
in the presence of such external conditions as temperature and chemical potentials, magnetic field, finite size effects 
etc (see, e.g., the papers @, H, 13: S S 0, B Hi and references therein). In particular, the color superconductivity 0,111 
as well as parity vaiolation and charged pion condensation EH El, El, EH EH] phenomena of dense quark matter 
, were investigated in the framework of these QCD-like effective models. 
t-H ' It is necessary to note that an effective description of QCD in terms of NJL models, i.e. through an employment 
of four-fermion theories in (3+l)-dimensional space-time, is usually valid only at rather low energies and densities. 
Besides, at present time there is the consensus that another class of theories, the set of (l+l)-dimensional Gross-Neveu 
' (GN) type models [H, E3], can a ^ so be used for a reasonable qualitative consideration of the QCD properties without 
C^) . any restrictions on the energy/density values, which is in an encouraging contrast with NJL models. Indeed, the 
GN type models are renormalizable, the asymptotic freedom and spontaneous chiral symmetry breaking are another 
properties inherent both for QCD and GN theories etc. In addition, the fi — T phase diagram is qualitatively the same 
in QCD and GN model [H, EH HO, HH H3] (here /i is the quark number chemical potential and T is the temperature). 
Note also that GN type models are suitable for the description of physics in quasi one-dimensional condensed matter 
systems like polyacetylene [23j . Thus, due to the relative simplicity of GN models in the leading order of the large 
7V c -expansion (7V C is the number of colored quarks), their usage is convenient for the application of nonperturbative 
methods in quantum field theory [241 ]. 

Before investigating different physical effects relevant to a real (3+l)-dimensional world in the framework of two- 
dimensional GN models, let us recall that there is a no-go theorem by Mermin-Wagner-Coleman forbidding the 
spontaneous breaking of continuous symmetries in two dimensions (25j . This theorem is based on the fact that in 
(l+l)-dimensional spacetime the Green function (correlator) of two scalar fields has at large distances a behavior 
| a; — yl -1 /^. Thus, if we take the limit \x — y\ — > oo first, the correlator vanishes at finite N c and, according to 
the cluster property, we formally obtain a zero vacuum expectation value of the scalar field, i.e. a prohibition of 
spontaneous symmetry breaking. However, there is a way to overcome this no-go theorem. Indeed, if the limit 
N c — > oo is taken first, then for \x — y\ — > oo we formally obtain a nonzero vacuum expectation value for the scalar 
field, i.e. the possibility for spontaneous symmetry breaking. It means that just the leading order of the large 
N c approximation supplies us in any (l+l)-dimensional model with a consistent field theory in which spontaneous 
symmetry breaking might occur. At present time this fact is well understood (see, e.g., the discussion in [20l l2ll |23| ) . 
This result restricts the range of validity of the no-go theorem to the finite iV c -case only. Clearly, since the no-go 
theorem does not work in the limit N c — > oo, the investigation of any low-dimensional model in the leading order at 
N c — > oo is much more physically appealing than the consideration of the model at finite N c . 

By this reason, such phenomena of dense QCD as color superconductivity (spontaneous breaking of the color sym- 
metry) or charged pion condensation (spontaneous breaking of the continuous isospin symmetry) might be simulated 
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in terms of simpler (l+l)-dimensional GN-type models in the leading order of the large N c approximation (see, e.g., 
PH and correspondingly). 

In our previous paper [26| the phase diagram of the (l+l)-dimensional GN model with two massless quark flavors 
was investigated under the constraint that quark matter occupies a finite space volume (see also the relevant papers 
[27|). In particular, the charged pion condensation phenomenon in cold quark matter with zero baryonic density, i.e. 
at fi = 0, but nonzero isotopic density, i.e. with nonzero isospin chemical potential /i/, was studied there in the large 
./Vc-limit. In contrast, in the present paper we consider, in the leading order of the l/iV c -expansion, the phase portrait 
of the above mentioned massive GN model in a more general case, where, for simplicity, temperature is taken to be 
zero, but both isospin and quark number chemical potentials are nonzero, i.e. fxi 7^ and fi ^ 0, and spacetime is 
considered to have the usual topology, R 1 x R 1 . Our consideration is based on the case of homogeneous condensates 
(an extension to inhomogeneous condensates in the case of /!/ = was recently considered in [22|, [HI). We suppose 
that these investigations will shed some new light on the physics of cold dense and isotopically asymmetric quark 
matter which might exist in compact stars, where baryon density is obviously nonzero (i.e. fj, ^ 0). 

The paper is organized as follows. In Sections II and III the effective action and thermodynamic potential of the 
two-flavored massive Gross-Neveu model are obtained in the presence of quark number as well as isotopic chemical 
potentials. In Section IV the phase structure of the model is investigated both in different particular cases (ft ^ 0, 
/ij = etc) and in the general case of [i ^ 0, fii 0. It turns out that at fij ^ and rather small values of /Li, the 
gapped pion condensed phase (PC) occurs. However, at larger values of /i several normal dense quark matter phases 
(without PC) are found to exist with different quasiparticle excitation properties of their ground states. In Section V 
the meson mass spectrum of each phase is discussed. Some technical details concerning the effective action and quark 
propagator are relegated to two Appendices. 



II. THE MODEL AND ITS EFFECTIVE ACTION 



We consider a (l+l)-dimensional model which describes dense quark matter with two massive quark flavors (u and 
d quarks). Its Lagrangian has the form 



L = q 



u-a i i I 1 ! 

7 \d v - to + /i7 + -y T37 



G 

q+ N, 



{qqf + (qrfrqf 



(1) 



where the quark field q(x) = qi a {x) is a flavor doublet (i = 1, 2 or i = u, d) and color V c -plet (a = 1, N c ) as well as 
a two-component Dirac spinor (the summation in ([T]) over flavor, color, and spinor indices is implied); Tk (k = 1, 2, 3) 
are Pauli matrices; the quark number chemical potential \l in ([1]) is responsible for the nonzero baryonic density of 
quark matter, whereas the isospin chemical potential /ij is taken into account in order to study properties of quark 
matter at nonzero isospin densities (in this case the densities of u and d quarks are different). Evidently, the model 
((T|) is a simple generalization of the original (l+l)-dimensional Gross-Neveu model [l6| with a single massless quark 
color V c -plet to the case of two massive quark flavors and additional chemical potentials. As a result, in the case 
under consideration we have a modified flavor symmetry group, which depends essentially on wether the bare quark 
mass mo and isospin chemical potential /ij take zero or nonzero values. Indeed, at jji = 0, too — the Lagrangian ([TJ) 
is invariant under transformations from the chiral SUl{2) x SUr(2) group. Then, at /ij ^ 0, too = this symmetry 
is reduced to Uj 3 l(1) x Ui 3 r(1), where I3 = T3/2 is the third component of the isospin operator (here and above the 
subscripts L,R mean that the corresponding group acts only on left, right handed spinors, respectively). Evidently, 
this symmetry can also be presented as Ui 3 (l) x Uai 3 (1), where Ui 3 (l) is the isospin subgroup and £7a/ 3 (1) is the 
axial isospin subgroup. Quarks are transformed under these subgroups as q — > exp(iar3)<7 and q — > exp(ia7 5 T3)g, 
respectively. In the case too 7^ 0, Hi = the Lagrangian ([TJ) is invariant with respect to the SUi(2), which is a diagonal 
subgroup of the chiral SUl{2) x SUr(2) group. Finally, in the most general case with too ^ 0,/!/ 7^ the initial 
model ([1]) is symmetric under the above mentioned isospin subgroup C// 3 (l). In addition, in all foregoing cases the 
model is color SU(V C ) invariant. 

The linearized version of the Lagrangian ([T]), which contains composite bosonic fields a(x) and Tr a (x) (a — 1,2,3), 
has the following form: 



L 



7 idy - to + ^7 + -yT37 



a - 17 TT a T a 



El 

4G 



00 + ir a ir 



From the Lagrangian @ one gets the following constraint equations for the bosonic fields 

G 



a{x) 



-2-^(<?g); 7r Q (x) 



G 



(2) 



(3) 



Obviously, the Lagrangian @ is equivalent to the Lagrangian ([I} when using the constraint equations Fur- 
thermore, it is clear that the bosonic fields ([3]) are transforming under the isospin Ui 3 (l) subgroup in the following 
manner: 



Ui 3 (l) : a — > cr; 773^-^3; 7Ti — > cos(2a)7Ti + sin(2a)7T2; 1^2 —> cos(2a)7T2 — sin(2a)7Ti, 



(4) 
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i.e the expression (tt J + w 2 ) remains unchanged under an action of the isospin subgroup Uj 3 (l). 

There is a common footing for obtaining both the thermodynamic potential and one-particle irreducible Green 
functions of bosonic <t(x) and ir a (x) fields ([3|) which is based on the effective action S c h (a, ir a ) of the model. In the 
leading order of the large iV c -expansion (corresponding to the one fermion-loop or mean field approximation), this 
quantity is defined in terms of the Lagrangian ([2]) through the relation 

exp(i5 cff (CT,7r a )) = N' J[dq][dq] exp(i y*id 2 ^, (5) 

where N' is a normalization constant. It is clear from and ([5]) that 

r 2 i 2 

S cS {v,ir a ) = -N c j ° 4 J a d 2 x + S oS , (6) 
where the quark contribution to the effective action, i.e. the term S e g in ([5]), is given by: 

exp(i5 cff ) = N' J [dq][dq] exp(i J qDqd 2 x^j = detD. (7) 

Here we used the notations 

D = \rfd v - m + m° + ^37° -v- rfn a T a (8) 

and v = m/2. Note also that D is a nontrivial operator in coordinate (x), spinor (s), and flavor (f) spaces, but 
it is proportional to the unit operator in the iV c -dimensional color (c) space. Then, using the general formula 
dot D = exp Tr xc f s In D, one obtains the following expression for the effective action: 

/2 _i_ 2 
a 4 J a d 2 x - iN c Ti sfx In D, (9) 

where we have taken into account that the trace of the operator lnD over the color space is proportional to N c . 
Starting from ([5]), one can define the thermodynamic potential (TDP) of the model in the mean-field approximation: 

S cS = -N c Q^(a,TT a ) ( d 2 x. (10) 

{(j.-Ha— const) J 

The ground state expectation values (mean values) of the bosonic fields, (<r(x)) = a° and (ir a (x)) = tt°, are solutions 
of the gap equations for the TDP Cl^^a, 7r a ) (in our approach all ground state expectation values do not depend on 
coordinates x): 

0fc = o ^1 = 0, where a = 1,2, 3. (11) 

OCT OTT a 

In particular, it follows from (jTTj) that if mo ^ then Trf = 0. In addition, one can put 7r| = 0, since the effective 
action depends on tk\ and W2 fields through the combination (irf + tt 2 ). Next, let us perform the following shift of 
bosonic fields in ([9]): a(x) — > <r(x) + <j°, iri(x) — > tt\{x) + n°, whereas the other bosonic fields, 7T2,3, stay unshifted. 
(Obviously, after shifting the new bosonic fields <j{x),ir a (x) now denote the small quantum fluctuations around the 
mean values a°, tt° of mesons rather than the original fields ([3J)). Moreover, we use the notations a° = M — mo and 
7r° = A. In this case 



D 



(vfdv — M + M7° + ^37° - i7 5 An) - (a(x) + i 7 5 ^(x)r a ) = 5,7 1 - (a{x) + i 7 5 ^ a (x)r a ) , (12) 



where So is the quark propagator which is a 2x 2 matrix in the flavor space, presented in Appendix|B] Then, expanding 
the obtained expression into a Taylor-series up to second order of small bosonic fluctuations a(x),ir a (x), we have 

S cS (a, 7r a ) = + S$ (a, tt ) + • • • , (13) 
where (due to the gap equations, the linear term in meson and diquark fields is absent in (TH 
1 

1 -rai , , f ,0 <7~ z + i 



N, 



(0) = _J d2x (M rn^ + A 2 ^ ^ _ %A ^ A) J ^ (w) 

S^(a,7r a ) = - [ d 2 x?-^ + ^Tr sfx {S (a + rfn a T a ) S Q (a + h 5 n a T a )} ■ (15) 
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The TDP fl^ tU (M, A) from (|14[) will be calculated in the next section, where on the basis of this function the phase 
structure of the GN model (1) in the leading order over 1/N C is considered. Note also that in (flU|) and (fT5)) the bosonic 
fluctuation fields a, ir a are really the coordinate dependent quantities. The trace of the So-operator and the products 
of er, 7r Q -fields in (fT5|) should be understood in the sense of formula (|A2[) (see Appendix [5} . Note the remarkable 
property that the effective action (fl~5|) is a generating functional of two-point and one-particle irreducible (1PI) Green 
functions of a- and 7r-mesons. Indeed: 

S 2 S (2) 

T X Y(x-y) = - xv , °* y (16) 

where X(x),Y(x) — a(x),ir a (x) and Txy( x — y) is the 1PI Green function of the fields X(x),Y(x). Variational 
derivatives in (TTB]) should be taken in accordance with the general formula (|A3[) (see Appendix [XJ. In the following, 
on the basis of these Green functions we study the meson mass spectrum in different phases of the model. 

III. THERMODYNAMIC POTENTIAL 

The Fourier transformation Sq 1 (p) of the inverse quark propagator S^ 1 (|12|) has the form: 

S^ip) =jrf + M7° + ^ 37 ° — M — i 7 5 An. (17) 
Clearly, in the direct product of spinor and flavor spaces it is a 4x4 matrix, which has four eigenvalues: 



ei )2 ,3,4 = -M ± V (po + /i) 2 - pj - A 2 + ^ ± 2zV(Po + ^) 2 -A 2 . (18) 
Then, applying the general formula (|A5|) to the expression p4[) for the thermodynamic potential, one gets: 
o ,ua ( (Af-m ) 2 + A 2 .A f d 2 p , . . 

jW(m,a) = ^ ^ + _Ji ln(ei) 

(M - m ) 2 + A 2 . /" d 2 p 



4G J (2tt) 



+ i/ -^ln{[(po + tf-(E+f][(p +rf-(E A f}}, (19) 



where = */ (i? ± ) 2 + A 2 , = E±v, z/ = /xj/2 and E = \/f>f + M 2 . The system of gap equations directly follows 
from (fTU): 



Q = dfl ll ^(M, A) _ M - m _ ^ /" d 2 p f E+ £?" 



g P f 

~ 2G J (2^) 2 El(p + M ) 2 -(E+) 2 ' ( Po +m) 2 -(^a) 

g^(M,A) A_9iA / ^ f 1 i 1 ] 

dA ~ 2G J {2*y\{p + rf-{E+y ( Po + n) 2 -(E A y)- 



(20) 



The TDP n^v^M, A) is symmetric under the transformations /i — > — fj, and/or /i/ — > — /ij. Hence, it is sufficient to 
consider only the region /i > 0, /// > 0. In this case, one can integrate in (|19p over po with the help of the formula 



^ln 
2tt 



(p + a) 2 -b 2 ] =L{\ a -b\ + \a + b\} 



(which is valid up to an infinite constant independent of quantities a, b) and obtain: 

\2 i A2 



n M ,.(M, A) = {M ^ +A - jf ^{e+ + E A + (/i - E+)9(n E+) + ( M - E A )9(v E A )}, (21) 



where 0(a;) is the Heaviside theta-function. In a similar way, the system of gap equations (|20|) is transformed to the 
following one: 

n _ dEV(M,A) _ M-mp f 00 d Pl f fl(E+-^)E+ fl(E A - M )E~ 1 

°~ " 2G M J_ 00 2^E\ E+ + E A /' 1 j 

n _ afi MlI/ (M, A) _ _A _ A /°° d Pl ; 0(E+ - M ) , 0(E A - M ) 



9A 2G ./ . 2„ . ,,, A ,. A 



ttI E+ E7 J 1 J 



The coordinates (gap values) M and A of the global minimum point of the TDP (|21[1 supply us with two ground 
state expectation values (qq) and (qi^Tiq), respectively, through the relations M = tuq + (a), A = (ni) and formulae 
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([3]). In particular, if the gap A is equal to zero, the ground state of the model is isotopically symmetric and there 
is no condensation of charged pions. However, if A ^ 0, then the ground state describes the phase with charged 
pion condensation, where the isospin Ui 3 (l) symmetry is spontaneously broken. In this phase the space parity is also 
spontaneously broken. Note also that the physical essence of the other gap M is the dynamical quark mass which is 
not equal to the bare mass mo, evidently. 

It is clear that the TDP (|21|) is an ultraviolet divergent quantity, so one should renormalize it, using a special 
dependence of the bare quantities such as the bare coupling constant G and the bare quark mass mo on the cutoff 
parameter A (A restricts the integration region in the divergent integrals, \p\\ < A). The renormalization procedure 
for the simplest massive GN model was already discussed in the literature, see, e.g., in [lj| [2(| H^] ■ In a similar way, 
it is easy to see that, cutting of the divergent integral in (|2"Tj) and using the substitution G = G(A) and m = mG(A), 
where 



1 



2G(A) 




(24) 



and m is a new free finite renormalization invariant massive parameter 1 (which does not depend on the cutoff A), it 
is possible to obtain for the TDP (|2ip a finite renormalization-invariant expression. Namely, 



n^ v {M,A) = lim ^(M,A;A), 



A^oo 



where 

&n,v(M, A; A) 



M 2 



mM 



4G(A) 



-A 



2tt I A 



E A + (/x - Efiefa - E+) + E A )0(n - E A )} 



A 2 

7T 



(25) 



(26) 



(To obtain (j2"B"|) we have omitted the unessential constant ^ as well as have added another one, ^-.) In (|2~4"|) the 
cutoff independent quantity Mq is the dynamically generated quark mass in the vacuum, i.e. at jj, = and /ir = 0, 
taken in the chiral limit, i.e. at mo = (see below). (The renormalized expressions for the gap equations are obtained 
in the limit A — > oo, if the replacements G — > G(A), mo — > mG(A) and \p\\ < A are done in (|23|) . or by a direct 
differentiation of the expression (|25j) .) The expression ([26| can also be presented in the alternative form 



n fl< „(M,A;A) = V {M,A;A) 



mM 



-A 



2tt l A 



E\ 



Pi 



M 2 



-(/x - EDe^ - E+) + ( M - E±)9(ji - E A )}, 



where 



F (M,A;A) 



M 2 + A 2 



4G(A) 

Obviously, the integral in (|27|) is convergent at A — * oo. Since 



dpi 



pi 



M 2 + A 2 + — 



A 2 

7T 



lim V (M,A;A) 

A^oo 



M' 2 



,2 r 



2tt 



In 



M 



f2 



M 2 



1 



V (M,A), 



(27) 



(28) 



(29) 



one can easily obtain from ([25]) . (|27p . and (|29]) the following finite renormalization-invariant expression for the TDP: 

mM 



Q litV {M,A) = V (M, A) 



J-oo 2tt I 



£1 



2^/p 2 



M 2 



-(/i - £;+)0( M - E+) + (fi— E A )0(p - E A )}. 



(30) 



Note that the integral in (|30|) is convergent. In the particular case of /x = 0, /i/ = and m = 0, i.e. for the massless 
GN model in the vacuum, it follows from (1301): 



n M , v (M,A) 



M 2 + A 2 


r /m 2 + a 2 \ 








- I 


fi=0,i/=0,ro=0 27T 


rl Mo 2 J 



(31) 



Note, the quantity m does not equal to the physical or dynamical quark mass M. The last one is defined by the pole position of the 
quark propagator. Alternatively, it can be found as a gap, i.e. one of the coordinates of the global minimum point of the thermodynamic 
potential (see also the remark in the paragraph just after H23H ), 
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Since for a strongly interacting system the space-parity in the vacuum is expected to be a conserved quantity we 
put A equal to zero in (|3"Tj) . As a result, the global minimum of the TDP (|3"Tj) lies in the point M = M , which 
means that in the vacuum and at mo = the dynamically generated quark mass is just the parameter Mo introduced 
in (|24[) . However, in the general case, i.e. at nonzero values of the chemical potentials, the dynamical quark mass 
depends certainly on /i, fij and obeys the system of the gap equations (|22|) - (f23]) (or (|20|) ). Another free parameter of 
the massive GN model, the quantity m, is not directly related to the quark mass, but rather to the mass of 7r-mesons. 

In the following, when studying the phase structure or the meson mass spectrum, the quantity Mq is still treated 
as a free parameter, however the massive parameter of the model, m = aMg/n, is fixed by a = ctQ ~ 0.17. In this 
case the vacuum properties of the massive GN model resemble the situation in some NJL-type models in realistic 
(3+l)-spacetime (for a more detailed discussion, see the next section ITV Gp . 

For the forthcoming investigations we need also the expressions for the density of quark number n q and isospin 
density nj: 



~ = f ^{^-E+) + 6(n-E A )}, (32) 

f-i££{^-'>-^-4 < 33 > 



IV. PHASE STRUCTURE OF THE MODEL 



A. Particular case: /i = 0, fii — 

Introducing the notation m = a.Ma/%, one can get from (|3"D"|) the following expression for the TDP at /i = 0, Hi = 
(usually, this quantity is called effective potential): 



Af 2 + A 2 
Q o( M,A) = ^t^- 

The corresponding gap equations look like: 



/M^+A 2 



Z7T 



*WA) 5!Mll f^|.^ = D| (35) 



DM ™ V M o 

2^o(M,A)^ 2Aln /M 2 4^X 

dA \ Af 2 / v ; 

The gap system (f35|) - ([36|) has several solutions, but the global minimum point (GMP) of the TDP (|34[> corresponds 
to the value A = 0. Then, at A = 0, the equation (|33|) vs M has three solutions of different signs. Just the one with 
largest absolute value corresponds to the GMP of the TDP. This quantity (gap) is denoted by M and depicted in 
Fig. 1 as a function of the variable a. Since the quark number density n q and isospin density m (|32p - (l33[) are equal 
to zero in this GMP, we conclude that at [i — and hi — the ground state of the model corresponds to the empty 
space, i.e. to the vacuum. Hence, in this case the gap M is the dynamical quark mass in the vacuum. Clearly, the 
gap M coincides with Mo in the chiral limit, a = 0. In addition, in Fig. 1 the behavior of the 7r-meson mass M w vs 
a in the case of fi = and hi = is also presented (it is the solution of the equation ([57)) from section IV A From 
the investigations of section IIV CI it will become clear that M^ coincides with the critical value /i/ c of the isotopical 
chemical potential /i/, at which the system passes from the vacuum state to the pion condensed phase. Just this fact 
is reflected in Fig. 1. Moreover, we have also depicted in this figure the behavior of the critical value fi c vs a of the 
chemical potential [i, at which the system passes from the vacuum to the normal quark matter phase at v = (see 
the next section ITV Bp . 

It is easily seen from Fig. 1 that the relation between the gap M in the vacuum and the pion mass M T (at fi = 
and [ii — 0) has a strong a-dependency and for some values of this parameter does not describe real physics. Recall, 
in real (3+l)-dimcnsional physical models the dynamical quark mass M is usually greater than M^ at fi = and 
Hi = and depends on the model parameters (coupling constants, cutoff parameter etc). In particular, the values 
M = 350 MeV and M^ — 140 MeV, i.e. M/M n = 5/2, are often used in the NJL-model investigations of dense quark 
matter [301 ] . So, in the following consideration of the phase structure of the model (1) and its meson mass spectrum 
in the most general case of fi ^ and v ^ 0, we will suppose the same relation between M and M n at fi — and 
Hi = 0. Evidently (see Fig. 1), this choice corresponds to a = ao w 0.17. Having fixed the parameter a — ao w 0.17, 
it is then possible to obtain M/Mq « 1.04, M^/Mq w 0.42 and m/Mo « 0.05, where Mq is the dynamical quark mass 
in the massless GN model at fi — and (j>i = 0. 



7 




FIG. 1: Dynamical quark mass M (curve 1) and 7r-meson 
mass Mir (curve 2) vs a = ■Km/Ma at fi — 0, fii = 0. The 
curve 3 is the critical value /j, c of the vacuum - normal 
quark matter phase transition (at fii = 0); the critical 
value Hie of the vacuum - PC phase transition is also 
given by the curve 2, i.e. fii c = M n (see in the section 

MM- 



FIG. 2: Dynamical quark mass M vs fi at /ii 
a = a » 0.17. Here fx c /Mo « 0.76. 



and 



Particular case: /i 0, fii — 



Using again the notation m = aMo/ir, one can get from ([30| the following expression for the TDP at fj, ^ 0, \ii = 0: 



n M (M, A) 



M 2 



aM M 



2tt 

6/(|U - VM 2 + A 2 ) 




/x^M 2 — M 2 — A 2 



(37) 



It follows from the gap equations for the TDP (|37p that A = in its global minimum point, whereas the M-coordinatc 
of the GMP obeys the equation: 



e^ 2 



M 2 )ln 



(p + V^ 2 - M 2 ) 2 
M 2 



6> (M 2 



Af 2 
M)ln M 



aM 
2M ' 



(38) 



Studying the GMP of the TDP (f3"T|) with the help of the stationary equation (|3"8")) . it is possible to show that at ju < /i c 
the GMP is arranged in the point (M, A = 0), where both the critical value /i c and the gap M are depicted in Fig. 
1. In this case the system is arranged in the vacuum state with n q = and ni = 0. However, if fi > /x c then the 
phase which is usually called the normal quark matter phase is realized in the model. In this phase the quark number 
density n q is nonzero, however the isospin density nj — at \xi = 0. In the particular case with a — «o ~ 0.17 the 
behavior of the M-coordinate (gap) of the GMP is presented in Fig. 2, where fi c w 0.76Mq, as a function of \i. 



C. Particular case: fi = 0, [ii 7^ 



In this case the TDP (|3"0")) has the following form: 

aM M 



Q V (M, A) = Vq(M, A) 



2tt 



dp 



-{ 



E A + E A 



2a/ P 2 + M 2 + A 2 }. 



The corresponding system of the gap equations looks like: 



27rdQ„(M,A) 
dM 

2irdQ u (M, A) 
9A 



2Mln 



2A 



hi 



M 2 + A 2 
M 2 

M 2 + A 



aM - 2M 



E- 



dpi\ 



EE+ EE A + Af 2 + A : 

1 2 

0, 



1^1 S A VpI + Af 2 + A 2 ■ 



(39) 

:}=0, (40) 
(41) 
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FIG. 3: The gaps A (curve 1) and M (curve 2) vs u = FIG. 4: Isospin density ni vs u at fi = and a = c*o ~ 

(j.i/2 in the case fi = and a = «o ~ 0.17. 0.17 in the vacuum (at v < v c « 0.21Mo), where m = 0, 

and in the PC phase (at f > v c ), where ni 7^ 0. 



where we have used the notations adopted after formula (fT^|) . The coordinates (gap values) M = M(v) and A = A(i/) 
of the global minimum point of the TDP (|39[) obey the gap equations (|40| and (|4Tj) . (In the present section we find 
it convenient to stress explicitly the fact that the GMP is indeed a function of the parameter v.) 

Recall the situation in (3+l)-dimensional NJL models with pion condensation, if the bare (current) quark mass is 
nonzero (l3j . In this case at some critical value /i/ c of the isospin chemical potential, which is just the pion meson 
mass M„ in the vacuum at /x = and /i/ = 0, i.e. /i/ c = M v , there is a continuous 2nd order phase transition from the 
vacuum phase (which is realized at v < v c = with M(v) = M(0) 7^ 0, A(^) = to the pion condensed one (at 

v > f c ), where M(v) 7^ 0, A(^) 7^ 0. This means that the TDP global minimum point (M(z/), A(z/)), corresponding 
to the pion condensed phase, has the following property: M(v) — » M(v c ) = M(0), A(^) — > 0, if 1^ — >■ v c+ . Here we 
again use the notations v = /ij/2 as well as M(0) for the dynamical quark mass in the vacuum. 

It turns out that the same qualitative picture of the pion condensed phase transition occurs in the framework of 
the massive GN model. Indeed, numerical investigations of the TDP ([55)) show that at some critical point v c there is 
a second order phase transition from the vacuum phase to the phase with charged pion condensation. It means that 
the GMP of the TDP is a continuous function vs v in the critical point v = v c . Now, in order to define v c and to 
prove that the equality v c — is also valid in the case of the massive GN model, it is necessary to remark that 

at v > v c the coordinates (M(u), A(i/)) of the GMP of the TDP (|39|) convert the expression in the square brackets 
of (|4Tj) into zero. Moreover, the equation (|40[) is also fulfilled. Since at v = v c we have a continuous phase transition, 
i.e. A(z/ C ) = 0, M(v c ) = M(0), 2 in the critical point v = v c this pair of equations is transformed into the following 



one 



M\Q) 



aM = 2M(0)ln— ^, (42) 

Mq 

T2( 



= 2^ / dpi ; , = , , — -• (43) 



M£ C J Q ^pi + M2(0)(pl+M2(0)- K 

Next, by inserting equation into the right hand side of (|4"2"j) , we find the useful relation 

aM _ 2 f°° , 1 



= 2^/ dPi ; 5 5 77^ ( 44 ) 



2M(0) c 7o v^TM2(0)(p 2 +Af 2 (0)-^ 

In the next sections we will study the meson masses in different phases of the model. In particular, we shall there 
derive an equation for the it- meson mass M n in the vacuum at fj, = and \ii = (comp. (|57p V Comparing (|44p with 
this equation, it follows that v c — i.e. the critical value \ii c is equal to the 7r-meson mass M v at /i = and 

fii = for arbitrary values of a. (Of course, one should take into account that the corresponding dynamical quark 
mass M appearing in this equation is nothing else than the parameter M(0) of the present section.) As a result, the 
dependence of p,i c and M n vs a is presented by the same curve 2 of Fig. 1. 



2 The quantity M(Q) vs a is nothing else than the gap M depicted in Fig. 1 as curve 1. 
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Clearly, at v < v c we have a phase which corresponds to the empty space (here both n q and nj are equal to zero). 
Due to this property, we use the notation vacuum for this phase. 3 In the vacuum phase one has A = 0, but the 
gap M is nonzero and does not depend on v (its behavior vs a is shown in Fig. 1). At v > v c the pion condensation 
(PC) phase with n q = and nj ^ is realized in the model. Inside this phase both gaps M and A are nonzero and 
depend on v. The isospin Ui 3 (l) symmetry is spontaneously broken in the PC phase. For the particular parameter 
value a — ao » 0.17 the behavior of gaps vs v is shown in Fig. 3, where v c as 0.21Mo. In Fig. 4 the isospin density 
nj vs v is presented. 



D. General case: fi / 0, /i; / 

In this case, starting from the TDP (|3H|) we obtain the following gap equations: 

dM \ A/ 2 1 



-2M 



r E+ E- 2 E+9{ji-EX) E~6(fi - E^) • 

I EE+ + EE A ~ ^/WT7& EE+ • 



2irdn^(M,A) = 
dA 



, M 2 + A 2 

111 1 ~w~ 



dpi \ V 1 , 2 g <"-5> g^-gj)\l =0 (46) 

£ A VS 2 + A 2 £ A J J 



Based on these equations, we have studied the properties of the GMP of the TDP (|30|) in the particular case of 
a = an w 0.17 and found the phase portrait, presented in Fig. 5. There the vacuum, pion condensation as well as 
three normal quark matter phases I, II and III are arranged. In the pion condensation phase the gaps A and M are 
nonzero, so here the isospin C/j 3 (l) symmetry is broken spontaneously. Throughout of this phase the gaps A and M 
do not depend on /i. It turns out that their dependencies on v in the PC phase at /i ^ are the same as in the PC 
phase at [i = (see Fig. 3). For points (y, fj), taken from the other phases of Fig. 5, the A-coordinate of the global 
minimum point of the TDP is zero (as a result, in these phases the isospin Ui 3 (l) symmetry remains intact), but the 
M-coordinate of GMP is not zero. Namely, inside the vacuum phase the gap M does not depend on {v,fj), i.e. it is 
a constant. (In particular, here M rs I.QAMq at a = ao ~ 0.17.) Our analysis shows that on the boundary between 
the vacuum and pion condensation phases the gaps are continuous functions vs (i and v. Hence, we conclude that a 
transition from the vacuum to the pion condensed phase or conversely is a second order one. 

It turns out that the gap M is a continuous (/i, ^-function inside each of domains I, II, and III of Fig. 5. In 
contrast, it is changed by a jump when each boundary between I, II, and III phases is crossed. To become convinced 
in this, look at Fig. 6, where the behavior of M vs v at two different fixed values of /x is presented (there, the phase II 
is shrank to the interval (a\, bi) at [i — 0.84Mo and to the interval (02, 62) at \i = 0.94Mq). As a result, we conclude 
that on these boundaries there is a first order phase transition. 

Now, let us consider the quark number density n q (|32p as well as the isospin density nj (I33[) inside each phase of 
the model. It is easy to see that for the vacuum phase these quantities are zero, thus justifying the name of these 
phase. Then, since the gaps A and M do not depend on fj, inside the pion condensed phase and the relations E^ > fi 
are true here, one can conclude that n q = in this phase and the isospin density ri/ vs v in the PC phase is presented 
in Fig. 4. For the normal quark matter phases I, II and III we have A = 0, so in order to obtain the expressions for 
n q and nj in these phases one can use the expression (13) of the paper [32[ for the quantity £l^ v {M, A = 0). As a 
result, in the phases I, II and III we have: 

n q = ^ + " - M) + v? - M 2 + ^-»\- Ml V { , _ v) 2 _ W rigaOi - (47) 
^ = g(g + v M) ^- - _ M2 _ 0(\fi -u\- M) ^ _ _ — 2 

TT IT 

Then, using the values of the gap M presented in Fig. 6, one can find the corresponding values of densities n q and 
nj shown in the curves of Fig. 7. It is clear from this figure that inside the Il-phase n q = nj. Since n q — n u + 



By definition, the vacuum is here the phase with zero densities n q and nj. However, one should realize that in a most general case the 
(dynamical) properties of its ground state depend on the values of fi and /ij. Indeed, in the model under consideration at /ij = there 
is an SUi(2) symmetry of the ground state in the vacuum phase. As a result, all three pions have a common mass. However, at /ij 7^ 0, 
i.e. when the ground state symmetry is reduced to the Ui 3 (l) subgroup, 7r-mesons have different masses in this phase (comp. section 
IVB2t . 
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FIG. 5: Phase portrait of the model in terms of /i and 
v = m/2 at a = a « 0.17. Here v c /M « 0.21, /x c /Mo w 
0.76. All lines of the figure are curves of first order phase 
transitions, except the boundary between the vacuum and 
PC phase. In the normal quark matter phases I, II and III 
pion condensation is absent (their properties are discussed 
in the text). 



FIG. 6: The behavior of the gap M vs v in the phases 
I, II and III at /j,/M = 0.84 (curve 1) and fi/M = 0.94 
(curve 2) at a = a « 0.17. Here ai « 0.47, a 2 ~ 0.62, 
6i « 1.11, b 2 ~ 1.19. 



and 2nj — n u — rid, where n u , rid are the densities of up and down quarks, correspondingly, it is clear from the above 
mentioned constraint that in the phase II the relation n u = —3rid is valid. 

Up to now we have studied thermodynamic properties of the model phases. Now the consideration of their dynamical 
peculiarities are in order. The first point we would like to discuss here is the spectrum of quasiparticles. In condensed 
matter physics they are simply the one-fermion excitations of the corresponding ground state. Recall, in the most 
general case the energy spectrum of u-, d-, u-, d-quasiparticles (quarks) are given in (|B6|) (see Appendix [Bj . It is clear 
from this formula that in the vacuum phase, where A = and M ~ 1.04Mo, the energy which is needed for a creation 
of the u- and d-quasiparticles is always greater than zero. Hence, both u- and c?-quarks are the gapped excitations 
of the vacuum phase. The similar property of a ground state is valid for the PC phase of the model, where also a 
finite amount of energy is needed to create up and/or down quarks. Due to this reason, the name gapped phases are 
usually used in these cases. 

However, in the case with normal quark matter phases I, II and III the situation is opposite. Indeed, it is easy to 
check that in the phase I both u- and d-quasiparticles are gapless. It means that there are no energy costs to create these 
quarks, i.e. there exist space momenta p* and p\* such that Pouip*) — and pod{p**) — 0, where pouipi) and Podipi) 
are the energies given in (|B6|) of corresponding quasiparticles. (For example, the point (y = 0.2Mq,[i = 0.84Mo) lies 
in the phase I with M w 0.088M , A = 0. Then it is easy to find from that p\ w 1.04M and p\* k, 0.63M o .) In 
contrast, in the phases II and III only u-quasiparticles are gapless, but (i-quarks are gapped. Note, some dynamical 
effects in dense matter such as transport phenomena (e.g., conductivities etc) depend essentially on the fact whether 
or not gapless excitations of the medium are possible. Hence, these effects can occur in a qualitatively different way 
in the phase I on one hand, and in the phases II and III, on the other hand. 

Finally, it is necessary to remark that the spectrum of mesonic excitations is also has a sharp phase dependence. In 
particular, in Fig. 8 the behavior of the 7r°-meson mass in the phase II is depicted at two different values of /i. It turns 
out that in this phase the 7r°-meson is a stable particle (at least with respect to strong interactions). However, in the 
neighboring phases I and III it is no longer a stable particle but a resonance. This fact as well as other peculiarities 
of the meson spectrum in different phases of the model is the subject of our consideration in the next section. 



V. MESON MASSES IN DIFFERENT PHASES 



As was noted in section [XT] (see the text after (fT5)) ), the effective action (fT5|) can be used for obtaining meson masses 
in different phases of the model. For this purpose, one should find from the outset all two-point 1PI Green functions 
(|16[) of meson fields. These 1PI Green functions are the matrix elements of the 4x4 meson matrix r(x — y). Then it 
is necessary to get the Fourier transformation T(p) of the meson matrix T(x — y) and find its determinant in the rest 
frame, where the two-component energy-momentum vector p has the form p = (po, 0). The equation 



dctr(po) = o 



(49) 
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FIG. 7: Quark number density n q (curve 1) and isospin 
density ni (curve 2) vs v at n/Mo = 0.84. Here ai « 0.47, 
bi « 1.11. 



FIG. 8: The behavior of the ir -meson mass M n o vs v in 
the phase II at n/M = 0.84 (curve 1) and fi/M = 0.94 
(curve 2). The values a» and 6j are given in Fig. 6. 



has in the plane of the variable p^ four (real- or complex-valued) solutions, one of them is the mass squared of the 
scalar cr-meson, whereas the other three solutions give the mass squared of the pseudoscalar 7r-mesons. 

Detailed investigations of the meson matrix T(po) show that its matrix elements of the form T^ 3 x(po) or Txtt 3 (po): 
where X = a, m , tt 2 , are equal to zero in all phases of the model, i.e. the matrix T(po) is a reducible one. This means 
that the neutral pseudoscalar meson, n° = tt^, does not mix with the other mesons, ir^ = (tti ± i7T2)/v2 or a. As a 
result, one root of the equation (|49|) can be found through the equation r W3ff3 (po) = 0, which supplies us with the 
mass squared of the 7r°-meson in different phases of the model. The other three meson masses are the zeros of the 
determinant of the reduced meson matrix, whose matrix elements are two point 1PI Green functions of the fields er, 
7Ti, and 7T2. 



A. The mass of n -meson 



The corresponding two-point 1PI Green function looks like: 

r^o(z) = - s ^£ = S M +iTTs { Sll (z) 1 5 S 11 (-z) 1 5 + S 2 2(zh 5 S 22 (-z) 1 5 

- S 12 (z) 1 5 S 21 (-z) 1 5 - S 21 (z) 1 5 S 12 (-z) 1 5 }, (50) 

where z — x — y, and the matrix elements Sij (z) of the quark propagator are presented in (|B5|) . Note, the expression 
((ST))) is valid for all phases of the model. Now, let us consider it in each phase. 



1. Vacuum and normal quark matter phase I: The case v = 0, fi > 



To illustrate the technique, which was elaborated in details in the framework of NJL-models with the color super- 
conductivity phenomenon [301 ] . we start from the most simple case of v = corresponding to the vacuum and the 
phase I only (see section HV Bl and Fig. 5). Since for these phases A = 0, the last two terms in ([5X)| . proportional to 
S\ 2 {z), vanish. The corresponding Fourier transformation of the expression (|50|) now looks like: 

r^o(p) = +iTr s y {Siib + £07 5 Sii(fc)7 5 + S 22 (p + k) 1 5 S 22 (k) 1 5 } , (51) 

where the Fourier transformations Sij (p) can be easily determined from (|B5j) . Using in (|51|) the rest frame system, 
where p = (po,0), and calculating the trace over spinor indices, we have 

where E = \Jk\ + M 2 , and the dynamical quark mass M is given by the value of the M-coordinate of the GMP of the 
thermodynamic potential. As was noted in Appendix [Bl in (|52p fco an d (fco + Pa) are correspondingly the shorthand 
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notations for fco + ie ■ sign(fco) and (ko + po) + is ■ sign(fco +Po), where e — > 0+. The /co-integration in ([52]) is performed 
along the real axis in the complex fc -plane. We will close this contour by an infinite arc in the upper half of the 
complex fco-plane. Taking into account the above-mentioned rule for the fco-integration, we have inside the obtained 
closed contour of the integral in ([52"]) four poles of the integrand which are located in the following points: 

(fco)i = —E — + is0(/x + E), (k ) 2 = E- (i + ie9(fx - E), 

(fco)s = -E-h-pq + je9(h + E), (k ) A = E-n-p + k6(n-E). (53) 

Since fi > and E > 0, the sum of the corresponding residues of the integrand function in these poles results in the 
following fci-integration in (|52p : 

To renormalize the expression (|5"i|) we use the gap equation (|2"2"|) at A = v = 0: 



where we took into account that ttiq/G = m. Substituting ([55]) into ([54]) and using the relation m = aM^/ir, we have 

Note that the quantity r„o T o(po) is a multi- valued function of the variable p\ which is analytic on some complex 
Riemann manifold described by several sheets. The expression in the right hand side of ([55]) defines r T o w o(p ) just on 
the first physical sheet only, which is the whole complex pp-plane, except for the cut p$ > AM 2 along the real axis. 

Recall that the mass squared M 2 of 7r°-mesons is the zero of this 1PI Green function vs p^. The zero should lie 
either on the real axis in the first sheet of the p^-plane (in this case it corresponds to a stable particle with real value of 
M^o such that < M?L < AM 2 ) or in the second sheet, corresponding to a resonance. Since at v — a mass splitting 
between n- mesons is absent (see also the remark in the footnote [3]) , throughout the section we use the notation M T 
both for the n°- as well as for the 7r ± -meson mass. 

It is clear from Fig. 1 (see also section HVB [I that in the vacuum phase at v = the relation /i < fj, c < M is valid 
for arbitrary a-values, so the theta-function in ([56]) is equal to unity. As a result, we see that in the vacuum the 
7r-meson mass satisfies the following equation: 

aM 9 f°° dki 



2M ^Jk\ + M 2 {Ak{ + AM 2 -Ml) 



Supposing that the quantity M in ([57]) is just the gap depicted in Fig. 1 as the curve 1, one can solve numerically 
this equation with respect to the variable M w . It turns out that the solution M n lies in the first sheet of the Riemann 
manifold and hence obeys the relation M 2 < AM 2 . The quantity M v vs a is shown in Fig. 1 as the curve 2. 

In contrast, in the case of the phase I at v = the corresponding 1PI Green function l[56]) does not has zeros in the 
first Riemann sheet of the variable i.e. there are no stable (at least with respect to strong interactions) 7r-mesonic 
excitations of the phase I ground state. In this phase all 7r-mesons are resonances. 



2. Vacuum and normal quark matter phases I, II and III: The case i/ / 0, /i / 

Technically this is a more complicated case, but the main ideas of calculations do not change. So, omitting technical 
details, one can obtain the following expression for the two-point 1PI Green function of 7r°-mesons in the rest frame: 

IW (po) = ^L- P 2 J™^± E{ J 2 _^ ) [B{E + v - M ) + sign(£; - u)6(\E v\ M )] . (58) 

It is also a multi-valued function of the variable p\ which is analytical on the same Riemann manifold, where the 
Green function ([56]) is defined. On the first Riemann sheet and at real values of p\ such that < p\ < AM 2 it looks 
like: 



- aM 2p pa pod(p, + v-M) Po\f{p 

1 o_o nn = — arctan — — H — — arctan — 

2ttM ^4M 2 - p 2 y/AM 2 - p 2 ^AM 2 - p 2 (/i + v) V '4M 2 - p 2 



P0 9(\,-u\-M) pp^-vr-M* 

V^M 2 ~ P 2 \ [i - v ym 2 - v 2 
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Let us, for example, again consider the case a — oto w 0.17, then M w 1.04Mo (see Fig. 2). In this case, for the /i 
and v values taken from the vacuum phase of Fig. 5, each theta-function in the expression (|58[) is equal to unity. As 
a result, the Green function (|58|) at v ^ coincides with the 7r° Green function ([56| at ^ = 0. Hence, in the vacuum 
phase the mass of the 7r°-meson does not depend on both v and [i. It takes the value M„o 0.42A/o hi the case of 
a = ao. In the general case of arbitrary a-values, the M^o-mass in the vacuum phase at v ^ is simply the pion 
mass at v — (see the line 2 of Fig. 1). 

One can easily check that the expression (I59|) turns into zero at some point of the interval < po < 2M only in the 
case when /j,) lies in the phase II (the corresponding value of po is the mass of the 7r°-meson). At some fixed values 
of n the behavior of M n o vs v is presented in Fig. 8 at a — oq. In contrast, in the phases I and III the expression (|59|) 
has no zeros in the interval < po < 2M. Hence, in these phases ir° is not a stable particle, but rather a resonance. 



3. The pion condensation phase 

Now, let us study the 7r°-mass in the PC phase, where both gaps A and M are nonzero. To obtain a compact 
expression for the two-point 1PI Green function r w o T o(po), it is again necessary to eliminate in (|50[) the coupling 
constant with the help of the gap equation (j23|) . i.e. to use the following relation 

1 f°° dki r 9{E+ - n) 9{E A - n) -i 



2G 2vr I E+ E A , 

Then, after tedious but straightforward calculations which are similar to that of section IV A 11 it is possible to find 

dki E£ + E* 



Clearly, the mass of ir° in the PC phase is equal to the isotopic chemical potential fij and does not depend on /i. 



B. The masses of a- and n -mesons 



As was noted above, to get the masses of a- and 7r ± -mesons, it is necessary to find the zeros (in the rest frame 
with p = (pojO)) of the determinant of the reduced meson matrix composed from two-point 1PI Green functions of 
these particles. Our calculations show that the Green functions are of the form T a7!l 2 (po) ~ A. So, in the vacuum as 
well as in the phases I, II and III there is no mixing between a- and TTi^-fields which leads to a further reduction of 
the meson matrix. Hence, to find the mass of the cr-meson in these phases, it is sufficient to investigate the separate 
equation T aa (po) — 0. The equation detn(po) = with 

V 1 7T2?T1 (PO) , 1 ,r 27 r 2 (PO ) J 

then supplies us with the masses of 7r ± -mesons. 



1. a -meson in vacuum and I, II, III phases 



In these phases A = 0. So, On the basis of the effective action (TTSl) and using the relation ([To]) and the methods of 
the previous section IV A 11 it is possible to obtain the most general expression for the two-point 1PI Green function 
of the (j- meson both in vacuum and in the I-, II-, III phases of the model 

F«r(Po) = ^ - (Po - 4M 2 ) ^± E ^ E l_ p2 J O(E + v M ) + sign(£ - v)6(\E -u\- M )] • (62) 

, where again E = \/k\ + M 2 . Let us now suppose that the pair of chemical potentials (fi, v) belongs to the vacuum 
phase of Fig. 5, where, evidently, M > fi+ v. In this particular case the expression in the square brackets of ([62]) is 
equal to 2, so 

-vac aM , 2 ,, f 2s f°° dk l 1 ,f.o\ 

2-kM Jo 7r E(AE z -pf ) ) 

It follows from (|63p that in the chiral limit, when a = and M ^ 0, the cr-meson is a stable particle with mass equal 
to 2M. However, at arbitrary small a > the zero of the Green function (f6"3")) . located at the point p\ = AM 2 of 
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the first Riemann sheet at a = 0, shifts to the second Riemann sheet, signalling thus that in the vacuum phase of 
the massive GN model the cr-meson is a resonance. It is quite reasonable that at small values of a the mass of this 
resonance is near 2M. 

Now remark that for values of /i and v from the regions I, II or III of Fig. 5 the square brackets of the integrand 
in ([62]) . cannot be negative. As a result, for all real values of p\ such that < p\ < AM 2 the Green function r CT(T (po) 
is a positive quantity, i.e. it cannot become zero. Thus, in the phases I, II and III of the model the cr-meson is also a 
resonance. 4 

2. 7r ± -mesons in vacuum and I, II, III phases 

The squared masses of 7r ± -mesons in these phases are given by the zeros of the equation detn(po) = in the 
PQ-plane, where II(po) is the matrix (|61[) . To find its matrix elements, it is convenient to use in the effective action 
(fl~5|) the new fields n ± (x) — (ni(x) ± \tt2{x))/V2 instead of the old ones, tti^x). Then, it is natural to define the 
corresponding Green functions 1^+^- (x — y) etc, where 

i2c(2) 

r + (x-v) = cS 

ott (y)dir + {x) 

etc. 5 The Fourier transformations of these Green functions are connected with the matrix elements of the matrix 
n(po) ([61]) by the relations 

iWibo) = rv 27r2 (po) = T [Ttt+tt- (po) + r w -,r+ (po)] , 

iW 2 (po) = -r W27ri (p ) = - [r 7r+7r - (po) - r^- T + (po)] . (64) 

Then, the determinant of the matrix ([6"1"[) looks like: 

detn(p ) =r w+w -(po) - Tn-n+ipo). (65) 
Our straightforward analytical calculations show that 

r^+TT- (po) = r ff o ff o(p + mj), r w -^.+(po) = r w o n .o(/i/-p ), (66) 

where T^OTrO is the 1PI Green function of the 7r°-meson, presented in (|58|) . Now suppose that at pg — M^o the Green 
function of the 7r°-meson turns into zero, when the chemical potentials (//, /i/) are fixed at some values in the vacuum 
phase or the I-, II-, III phases. Then, on the basis of the relations (|66f it is clear that T w + W - (po) ~ ((po + M/) 2 — -^»o) 
and r 7r - 7r +(p ) ~ ((m/ — Po) 2 — M^o)- As a result, we see that 

dctn(p ) ~ [(po + Ai,) 2 - M 2 ] • [(/i/ - po) 2 - M^o] = [p^ _ ( Mjr0 - vi) 2 } ■ [Po - (M v o + f ij) 2 }. (67) 

Hence, the zeros of the determinant ([6"T]) . i.e. the quantities M 2 + = (M n o — /i/) 2 and M 2 _ = (M n o + Hi) 2 , can be 
identified with the mass squared of 7r ± -mesons. 

3. a- and -mesons in the pion condensation phase 

As noted at the beginning of the present section, there arises a mixing between a and 711,2 fields in the PC phase of 
the massive GN model. Thus, to define the mesonic mass spectrum one should find all the zeros of the determinant 
of the meson matrix, composed of corresponding two-point 1PI Green functions of the form (|16[) . We have found an 
exact analytical expressions for these Green functions and have shown that the determinant has a zero in the point 
Pq = 0. (In order not to overload the paper with rather cumbersome formulae, we do not present here the expressions 
for these Green functions.) It means that in the PC phase there is a massless bosonic excitation. It can be treated 
as a Goldstone boson which is a consequence of the spontaneous breaking of the isospin £// 3 (l) symmetry in the PC 
phase. 



4 Strictly speaking, the found resonance character of a is here associated to the existence of <yij-thresholds enabling the meson decay into a 
(non-observable) qg-pair. Clearly, in order to model e.g. the confinement properties of "more realistic" two-dimensional QCD ,33), one 
should consider a more sophisticated GN-model incorporating some suitable prescription for quark confinement (see, e.g., [341). Within 
such a model one could then treat the cr-decay into observable pions, a — > tttt, which is, however, outside the scope of this paper. 

5 In the phases with zero gap A the Green functions of the form r + + (x — y) and T^- Jr _ (x — y) vanish. 
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It turns out that further information about mesons in the PC phase can be found in the chiral limit, i.e. at tuq — 0. 
Indeed, in this case the Green functions of the form T r77ri2 (po) are identically equal to zero, so that the er-meson 
does not mix with 7ri.2-fields. Moreover, it is possible to show that in the massless GN model the Green function 
T aa (po) coincides in the PC phase with the Green function r w o T o(po) (see (JBOJ))- Due to this relation we conclude 
that M a = M n o — [ii in the PC phase of the massless GN model. 

VI. SUMMARY AND CONCLUSIONS 

Recent investigations of the phase diagram of isotopically asymmetric dense quark matter in terms of NJL models 
show that their pion condensation content is not yet fully understood. Indeed, the number of the charged pion 
condensation phases of the phase diagram depends strictly on the parameter set of the NJL model. It means that 
for different values of the coupling constant, cutoff parameter, bare quark mass etc just the same NJL model predicts 
different numbers of pion condensation phases of quark matter both with or without an electric neutrality constraint 
(see, e.g., [H, EH). Thus, to obtain a more objective information about the pion condensation phenomenon of dense 
quark matter, it is important to invoke alternative approaches. One of them, which qualitatively quite successfully 
imitates some of the QCD properties (see also the Introduction), is based on the consideration of this phenomenon in 
the framework of asymptotically free (l+l)-dimensional GN models in the leading order of the large 7V c -technique. 

In the present paper we have studied the phase structure of the massive GN model (1) in terms of quark number 
(//)- as well as isospin {pi) chemical potentials in the limit N c — ► oo (for simplicity, the temperature has been taken 
to be zero). After renormalization (comp. section ITTT]) . this model contains two free parameters: Mo-the dynamical 
quark mass in the vacuum of the corresponding massless GN model and the renormalization-invariant quark mass 
m = rno/G = aMo/ir (see also the remark in footnote [1]) . In our considerations we often put a — «o ~ 0.17 in 
order to have the same relation between the dynamical quark mass M and the 7r-meson mass M v in vacuum, i.e. 
M/Mt; — 5/2, as used in some other NJL model parameterizations [3(J. Just at a — ao the phase portrait of the 
model is presented in Fig. 5 in terms of fi and v — fJ-i/2. 

Firstly, we have found that at T = the charged pion condensation phase of the GN model is realized inside the 
(noncompact) chemical potential region [ii > M^o, where [i is not greater than Mq/V^ and M^o is the vacuum mass 
of the 7r°-meson. In this phase the isospin Ui 3 (l) symmetry is spontaneously broken down and a massless Goldstone 
bosonic excitation of the ground state appears. Moreover, we have shown that the mass of the 7r°-meson in the PC 
phase is equal to the isospin chemical potential \ii . All one-quark excitations are found to be gapped particles in this 
phase. As a result, the quark number density n q is equal to zero in the PC phase. 6 The same properties of the PC 
phase is predicted in the framework of some NJL model parameterizations (see, e.g., in [l2l. IT3|). In contrast, in the 
NJL phase diagram the pion condensation phases occupy a compact region and for some parametrization schemes the 
gapless pion condensation might occur [ll|, [l2|, [HI [H| ■ 

Secondly, at rather large values of the quark number chemical potential fi we have found a rather rich variety of 
normal quark matter phases I, II, and III (see Fig. 5), in which the quark number density n q does not vanish (see Fig. 
7). In particular, it turns out that in phase I both u- and ci-quarks are gapless quasiparticlcs. Contrary, in phases 
II and III only it-quarks are gapless, whereas d-quarks are gapped. By this reason, dynamical effects in transport 
phenomena for dense quark matter (e.g., conductivities etc) can occur in a qualitatively different way in the phases I 
and II, III. We have studied also the 7r-meson mass spectrum of these phases and found that in the phase I and III 
the 7r-mesons are resonances. However, the phase II is the so-called "stability island" for 7r-mesons. Indeed, as it was 
shown by our numerical calculations, the 7r°-meson is a stable excitation of the ground state of this phase. Its mass 
vs v is depicted in Fig. 8. The 7r ± -mesons are also stable in this phase, but their masses are M„± = |M w o =p (see 
section IVB2P . (The same relation between tt°- and 7r ± -meson masses is also valid inside the vacuum phase of Fig. 

In conclusion, by using the above rather simple approach to the GN phase diagram, we have found a variety of 
phases with rather rich dynamical contents. A related interesting issue could be the extension of these investigation 
to inhomogeneous condensates [Hj]. We hope that our investigation of the phase diagram of the massive GN model 
will shed some new light on the phase structure of QCD at nonzero baryonic and isotopic densities. Obviously, a more 
realistic imitation of the QCD phase diagram requires to include also a nonzero temperature as well as a suitable 
confinement prescription for quark propagators [34j |. 



6 In the gapped phases, PC or vacuum phases, the relations > fj, are valid. Then, using [I: >2 l i . it is clear that n q = in these phases. 
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APPENDIX A: TRACES OF OPERATORS AND THEIR PRODUCTS 

Let A, B, ... are some operators in the Hilbert space H of functions f(x) depending on two real variables, x = (x°, x 1 ) 
In the coordinate representation their matrix elements are A(x,y), B(x,y), correspondingly, so that 



(Af)(x) = J d 2 yA(x,y)f(y), (A ■ B){x,y) = J d 2 zA(x, z)B(z,y), etc 

By definition, 

Tri = J d 2 xA(x,x), Ty(A ■ B) = J d 2 xd 2 yA(x,y)B(y,x), etc. (Al) 

Each function f(x) G H can be considered as an operator /, acting in this space, with matrix elements f(x)5(x — y), 
where 6 (x — y) is the two-dimensional Dirac delta-function. As a result, one can formally consider the trace of 
functions, their products as well as the traces of more complicated expressions, such as the products of operators and 
functions. Indeed, using the definition (|Aip we have 

Tr/ = Tr/ = 6(0) J d 2 xf(x); Tr(A/ 2 ) = Tr(A ■ / 2 ) - J d 2 xd 2 yf 1 (x)6(x - y)f 2 (y)S(y - x) 
= 5(0) / d 2 xf 1 (x)f 2 (x); Tr{i/} = Tr{i ■ /} = / d 2 xd 2 y A(x, y)f(y)6(y - x) = ( d 2 xA(x,x)f(x); 



Tv{AhBf 2 } = Tr{i ■ fr ■ B ■ f 2 } = J d 2 xd 2 yd 2 vd 2 uA(x 1 v)f 1 (v)S(v - y)B(y,u)f 2 (u)S(u - x) 

= I d 2 ud 2 vA(u,v)f 1 {v)B{v,u)f 2 (u). (A2) 



In particular, it follows from (|A2|) that 

W"^-* -WCri- (A3) 

of{x) dfi{x)Sf 2 { y ) 

Now suppose that A(x,y) = A(x— y), B(x,y) = B(x — y), i.e. that A,B are translationally invariant operators. Then 
introducing the Fourier transformations of their matrix elements, i.e. 

A(p) = J d 2 zA{z)e^ z , A{z) - J ^A{p)e-^, etc, (A4) 

where z = x — y, it is possible to obtain from the above formulae 

Tri = A(0) j d 2 x = J 0fA(p) J d 2 x. (A5) 

If there is an operator function F(A), where A is a translationally invariant operator, then in the coordinate representa- 
tion its matrix elements depend on the difference (x — y). Obviously, it is possible to define the Fourier transformations 
F(A)(p) of its matrix elements, and the following relations are valid (A(p) is the Fourier transformation for the matrix 
element A(x — y)): 

F[A){p) = F(A(p)); TrF(A) = J ^L F {A(p)) J d 2 x. (A6) 

Finally, suppose that A is an operator in some internal n-dimensional vector space, in addition. Evidently, the same 
is valid for the Fourier transformation A(p) which is now some n x n matrix. Let Xi(p) are eigenvalues of the nxn 
matrix A(p), where z = 1,2, .., n. Then 

TrF(l) = J ^LtrF(Z(p)) jd 2 x = Y,j ^^(p)) J d 2 x. (A7) 

In this formula we use the notation tr for the trace of any operator in the internal n-dimensional vector space only, 
whereas the symbol Tr means the trace of an operator both in the coordinate and internal spaces. 
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APPENDIX B: QUARK PROPAGATOR 



It is clear from (|12[) that the quark propagator So is the following 2x2 matrix in the two-dimensional flavor space: 

c _ ( Su > S\2 \ / D + , D12 \ mi , 

So= [s 21 , s 22 )-{d 21 , d^) > ( B1 ) 

where (the summation over a = 0, 1 is implied) 

£>± = i 7 Q 5 Q - M + (// ± i/)t°, D 12 = D 21 = -i 7 5 A. (B2) 
The connection between Sij and Dtj is the following: 

S u = [D+ - D 12 DZ 1 D 21 ] ^ , S 21 = -DZ 1 D 21 S lx , 

S 22 = [D- - D 21 D^D 12 ] -1 , Su = -D+ l D 12 S 22 . (B3) 
It is easy to establish the following relations: 

D ± = J e-*<"-» 1 { (po + p. - E T b°A + + (po + P + E±)7°A- } , 

( flir . = /^ e ,*-..{^_ + ^/l (B4) 
7 (2tt) 2 [ p + fx - ET PQ + + E ± J 

where E ± = E±v, E = ^p? + M 2 , A± = |(1± 7 ° (7l g +M) ), and 71 = -7 1 . 7 Using the relations flB3]), it is possible 
to obtain from (lB2l) and (lEil): 



_^P_ e -i P (x-») ( Po + M + £ QT Po + M ~ E+ 



■ - ' (2-) 2 \(po+^-(Exr ' A - + (po+tf-mr 7 ( 



(2^) 2 I (po + A 1 ) 2 - (^a) 2 (Po + M) 2 - (EX) 

621 ~ iA ,W e ibo+^) 2 -(^) 2 + ( P o+M) 2 -(^ A ) 2 r (B5J 

where A± - |(1± 7 ° (71 ^~ M) ) and Po in the integrand is a shorthand notation for po-H£-sign(p ), where e — > 0+. This 
prescription for the quantity po correctly implements the role of the quantities \i and [ii as the chemical potentials 
and preserves the causality of the theory [3l|. It is worth also to note the following useful relations: 

7 5 A ±7 5 = A ± , 7°A ± 7 = A T . 

The poles of the matrix elements (|B5[) of the quark propagator in the energy-momentum space give the dispersion 
lows for quasiparticles, i.e. the momentum dependence of the quark (pou, Pod) and antiquark (poa, Pod) energies, in a 
medium 

Pou = E A -fi, PQd = E^-fi, p oa = -(-^a +/■*)' Pad = -( e a (B6) 



Strictly speaking, the quantities po«, pod from (|B6j) are the energies necessary for the creation of quarks with momen- 
tum pi, whereas pou, Pod i s the energy necessary for the annihilation of antiquarks. 
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